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Domain Decomposition Methods for Space Fractional Partial 

Differential Equations* 

Yingjun Jiang'^and Xuejun Xu^ 


Abstract 

In this paper, a two-level additive Schwarz preconditioner is proposed for solving the al¬ 
gebraic systems resulting from the finite element approximations of space fractional partial 
differential equations (SFPDEs). It is shown that the condition number of the preconditioned 
system is bounded by C(1 -I- H/S), where H is the maximum diameter of subdomains and 5 is 
the overlap size among the subdomains. Numerical results are given to support our theoretical 
findings. 

Keywords, fractional differential equations, overlapping domain decompositions, preconditioners 

1 Introduction 

Space fractional partial differential equations have been wildly used to describe the supper- 
diffusion processes in the natural world (see m)- Let U denote a polyhedral domain in and 
AI{z) denotes a probability density function on where = {z € || 2:||2 = 1} and 11-112 

denotes the standard Euclidean norm. In this paper, we consider the following multi-dimensional 
SFPDE( [I5]) 

— f D‘1^u{x) M{z)dz + cuix) = f{x), X G U, (1.1) 

Jsd-^ 

where l/2<a<l,c>0 and which will be given later, denotes the directional derivative 
of order 2a in the direction z. We assume M is symmetric about origin, i.e., M{z) = M{z') if 
z, z' G satisfy z -h z' = 0, which means that the above SFPDE is symmetric. 

Actually, the equation (jl.ip is an appropriate extension from one dimensional problem 

- (p+ q xD^)u + cu = f, (1.2) 

and its corresponding developing equation can be used to describe a general super-diffusion process 
(see m for details), where denote Riemann-Liouville fractional derivatives. One 
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(1.3) 


special case of is 

d 

- -oo^x“ + <li XiD‘^^)u + CU = f 

i=l 

with Pi,qi € M satisfying pi = qi and pi + qi = 1. We may find that equation (ll.3|) can be obtained 
from a by taking 

d 

M = '^Pi5{z - Ci) + qi5{z + ei), (1.4) 

i=\ 

where is the zth column of identity matrix in and 6 denotes the Dirac function on 

Extensive numerical methods have already been developed for SFPDEs, like finite difference 
methods (see e.g., [aimiisiEsiES]), finite element methods (see e.g., [51M fMl[22] l. and spectral 
methods m- Due to the nonlocal properties of fractional differential operators, the most impor¬ 
tant issue for numerical computation of SEPDEs is how to reduce the computation costs. Some 
methods for the reduction have already been designed, like alternating-direction implicit methods 
(ADI) [l8l[30l[3T] , special iterative methods and multigrid methods [201136] . 

The discrete systems Ax = 6 of SFPDEs usually have the following characteristics: 1). the 
condition number of A increases fast, as the mesh becomes fine; 2 ). the coefficient matrix A is 
dense. As for iterative methods, two issues need to be concerned for efficiency: one is how to 
construct good preconditioners for the discrete system Ax = b, which may help us to save the 
iterative steps; the other is how to reduce the computation cost of each iterative step, for which 
some papers employ the multiplication of Toeplitz matrices and vectors with relog(n) computation 
complexity (see e.g., [T^[^[^[^TH54] 1. If the iterative step can be carried out in parallel, the 
efficiency of solving SFPDEs may be significantly improved. The parallelizable algorithms have 
been wildly used in numerical solutions for PDFs (see e.g. [2H]). Whereas, to the best of our 
knowledge, no parallelizable algorithms have been designed for SFPDEs. 

In this paper, we shall construct a two-level additive preconditioner for the discrete system 
Ax = b resulting from the finite element approximation of (ini), and then nse the preconditioned 
conjugate method (PCG) to solve it. The preconditioner we construct is almost optimal, i.e., 
the condition number of the preconditioned system is bounded by C(I + H/S), where H is the 
maximum diameter of subdomains and 6 is the overlap size among the snbdomains. Moreover, 
the preconditioner may be employed in parallel or each step of the PCG can be carried out in 
parallel. As a result, the whole numerical solution processes, including the generation of A, b and 
the multiplication of matrices and vectors, may be conducted in parallel. 

Without loss of generality, we focus on the case d = 2, namely, we consider the problem 
dEB) in R^. For A C R^, denote T^(A) the space of all measurable function u on A satisfying 
fj^(v(x))^dx < oo, and C^(A) the space of infinitely differentiable functions with compact snpport 
in A. Set 

{v,w)a = jj vwdxdy, ||?;||a = (u,i;)^^, 

and they are abbreviated as {v,w) and ||u|| respectively if A = R^. 

To simplify our statement, we make a convention here: function v defined on a domain A C R^ 
also denotes its extension on R^ which extends v by zero outside A. The constant C with or without 
snbscript shall denote a generic positive constant which may take on different values in different 
places. These constants shall always be independent of mesh sizes and numbers of subdomains. 
Following [35], we also use symbols <, > and in this paper. That ai < bi, a 2 ^ 62 and 03 63 
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mean that ai < Cibi, a 2 > (^262 and < 03 < <^363 for some positive constants Ci, 6 * 2 , C 3 and 
^^3- 

The rest of the paper is organized as follows: in section 2, the variational problem of dni) and 
its finite element discretization are described; in section [31 the two-level additive preconditioner 
for the SFPDEs is presented; in section jH it is proved that the preconditioner is almost optimal; 
finally, in section 5, the numerical results shall be given to support our theoretical findings. 


2 The model problem and its discretization 

In this section, we shall describe the SFPDEs in details, and then introduce its variational formu¬ 
lation and finite element discretization. 


2.1 The model problem 

Definition 2.1. [ 8 ] Let /r > 0, 0 € M. The /rth order fractional integral in the direction 2 : = 
(cos 0 , sin 0 ) is defined by 


roo 

D~^v{x,y) -.= Dg^v{x,y) = j 


r cos 9,y — T sin 6)dT, 


where T is the Gamma function. 


Definition 2.2. [ 8 ] Let n be a positive integer, and 0 € M. The nth order derivative in the 
direction of z = (cos 0 , sin 0 ) is given by 


Dgv{x,y) 


„ d . 5 

cos 0 ——h sm 0 — 
ox ay 


n 


v{x,y). 


Definition 2.3. [ 8 ] Let /r > 0, 0 G M. Let n be the integer such that n — 1 < ^ < n, and define 
a = n — fj,. Then the /rth order directional derivative in the direction of 2 ; = (cos 0, sin 0) is defined 
by 

D^v{x,y) := D^v{x,y) = D^D"%(x,y). 

If V is viewed as a function in x, Dq, are just the left and the right Riemamm-Liouville 
derivatives (see Appendix). 

Definition 2.4. [ 8 ] Assume that u : ^ M, ^ > 0. The /rth order fractional derivative with 

respect to the measure M is defined as 


D^^v{x,y) := D’^v{x,y)M{9)d9, 

where 5^ = [0-|-n, 27r-|-n) with a suitable scalar n, and M{9), which satisfies M{9)d9 = 1, is 

a periodic function with period 27r. Usually we take u = 0, if it causes no unreasonable expression 
(see (II.4p ). 

For u : ^ M, define differential operator La in as 

LaU = —D^^u + cu. 
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Denote D a polygonal domain in set 1/2 < a < 1. The model problem of this paper is to find 
tt : D —>■ M such that 

f LaU = f, in D, 

\ u = 0, on di}, 

where / is a source term and we assume that M{9) satisfies M{6) = M{9 + vr) for 0 G R, i.e., 
(12.ip is a symmetric problem. Here, we recall the convection made in Section [1) u also denotes 
its extension by zero outside D. 


2.2 The variational formulation and finite element discretization 

Definition 2.5. [27] Let ^ > 0, Tv{^i,^ 2 ) be the Fourier transform of v{x,y), |^| = Vsi+^- 
Define norm 

(l + |^|y/2|^„|| . 

Let := {u € L^(R^); < oo}. 

For V G FTq (D), we also denote ||u||j:^h(r 2 ) by ||'y||HM(n). It is known that iL^(R^) is a Hilbert 
space equipped with the inner product = ((1 + T'rc) andC^PM^) is dense in 

iL^(R^) (see m)- Since we employ the finite element discretization, the weak fractional directional 
derivative need to be introduced. Let denote the set of locally integrable functions on 

R2. 




Definition 2.6. [9] Given /r > 0, 0 G R, let u G L^(R^). If there is a function G such 

that 

Dg+^w) = {v^,w), Vtc G C'^(R^), 

then is called the weak /rth order derivative in the direction of 9 for v, denoted by D^v, i.e., 
Vf, = D^v. 

The weak derivative DgV is unique if it exists and the weak derivative coincides with the 
correspondent derivative defined in Definition 12.31 if v G (^“(R^). In the following, we use DgV to 
denote the weak derivative. 

Lemma 2.7. Let y > 0. For any v G 0 < s < y, and 0 G R, the weak derivative 

DgV exists and satisfies 

J^DgviCi,C2) = {27ri^icos9 + 27ri^2sm9yTv{^i,^2), ( 2 . 2 ) 

IlDgull < C'||u||j^m(k2 ). (2.3) 

Define the bilinear form B : x —>■ R as 

r2'K 

B{u,v) ■.= - / {DQU,DQj^^v)M{9)d9 + c{u,v). 

Jo 

Because M{9) = M(0 + 7r) for 0 G R, it is easy to check that B{v,w) is a symmetric bilinear form, 
i.e., B{v,w) = B{w,v) for v,w G Hq{Q). The variational formulation of (|2.ip (see mi) is to find 
u G Hq{LI) such that 

B{u,v) = if,v), yvGH^in). (2.4) 
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Here M is taken such that (12.41) admits a unique solution in Hq{Q,) (for the details, we refer 

to [8l[9]). 

We construct a quasi-uniform triangulation Fj:/ = of H with 

the diameter of ~ 0{H). 

Divide each Dj into smaller simplices Tj of diameter 0{h), such that F/j = {tj} form a finer 
triangulation of D. Denote Vh and Vh piecewise linear finite element function spaces defined on 
the triangulations Fj;/ and F/^ respectively. It is known that I 4 C Hq{0) C Hq{Q). 

The finite element approximation for (12.4p (the details please see [9]) is to find Uh € Vh such 
that 

B(uh,v) = (f,v), Vv € Vh, (2.5) 

where B(v,w) = — {DgV, Dg^^w)M{6)d9+c{v,w), M{9) is equal to a discrete form Pfc(I(0— 
9h) such that B{-, ■) is a symmetric bilinear form, and 

B{v,v) H(u,'w;) < ||u||H«(o)|k||H“(n), v,w € (2.6) 

Meanwhile 

/*27r 

/ M{9)d9 < 1. 

Jo 

Remark 2.8. The direct finite element discretization of (12.41) is 

B{uh,v) = {f,v), v€Vh. (2.7) 

But the discretization is hardly computed. So we use (ESI) instead of (EH, where B{-,-) is 
understood as the approximation to B{-, •). For the details, please refer to [9]. 

3 A two-level additive Schwarz preconditioner 

Take fh G 14 such that {fh,v) = if,v), Vv G I 4 and define a linear operator A : Vh ^ Vh 
satisfying 

{Av,w) = B{v,w), Vv,w &Vh- (3-1) 

Since B{v,w) is a symmetric bilinear form, by (|2.6p . we know that H : I 4 ^ 14 is symmetric 
positive definite with respect to (•, •), i.e., 

[Av, w) = {v, Aw), v,w & Vh', [Av, v) > 0 , 0 7^ u G I4. 


Then bilinear form 

{v,w)a ■■= {Av,w) 

also induces an inner product on I 4 . Set norm 

\\v\\a = {Av,v)^/'^, veVh. 

By ()2.6I) . we have 

II^IU ~ IbllH-(Q), Vu G 14 . 


(3.2) 
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The problem (|2.5I) can be restated as to find such that 

Auh = fh- 


(3.3) 


For the above equation, we shall construct our two-level Schwarz preconditioner and then use 
PCG method to solve it. 

Our preconditioner is designed by making use of the following overlapping domain decompo¬ 
sition O = where the subdomain O' contains coarse subdomain Oj, and satisfies 

the diameter of O' 0{H). 

Meanwhile the boundary of O' align with the mesh of triangulation F/i, and the distance from 
90( n O to Oj is greater than 6 , which is a positive constant measuring the overlapping size among 
the subdomians. Define subspaces V) of Vh as 

Vi = {v e Vh\v{x) = 0, X G O \ O'}, f = 1, 2,... , J, 

and let Vq = Vh- 

For our analysis, we regroup the subregions in terms of the coloring strategy (see e.g. [28]). 
By a minimal or good coloring, we group the subregions {0(} into classes, each of which has 
some disjoint subregions and can be regarded as one subregion. Exactly, decompose the index set 
{ 1 , 2 ,... , J} = with li satisfying that 0) n 0'^ = 0 for any /, A: € Ii{l 7 ^ A); for f = 1 ,..., 

define new subregions Dj = Ujg/.D'- and new subspaces Vi = 0 jg 7 . Vj. 

For each A € {0,1, 2,... , J}, we define some projectors Qk, Pk '■ Vh ^ 14 by 

{QkV, w) = (v, w), (PkV, w)a = {v, w)a, Vtc €Vk,v G Vh, 

and define the linear operator Ak '■ Vk ^ 14 by 

{AkV,w) = {Av,w), v,w GVk- 


It is not hard to verify that 

AkPk = QkA, k = 0,1,... , J. (3.4) 

To help our analysis, we define some projectors Pi : Vh ^ V, i = 1,..., , by 

{PiV, w)a = {v, w)a, Vrc GVi,v G Vh- 

Remark 3.1. Different from the integer order PDEs, it is interesting to see that Pi 7 ^ Yljeii ^j- 
Now, we are ready to present our two-level additive Schwarz pre-conditioner, i.e., 


.7 

Bh = '^Ar^Qi. 


i=0 


By (|3.4p . we have 


BhA = A-^QiA = Y A-^AiPi = Y P- 


i=0 


i=0 


i=0 


Define Ph '■= then the preconditioned system is 

PhUh = Bhfh- 


(3.5) 


In the next section, we shall prove the condition number of Ph is bounded by C(1 -|- H/6), 
where the constant C is independent of mesh size and the numbers of subdomains, but dependent 
of J^. 
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4 Condition number estimate 


We first introduce two interpolation norms and relevant Sobolev spaces (see e.g., 127]). Let A be 
a domain in For integer m, denote by || • the Sobolev norm of integer order m, i.e., 





1/2 


D^v 


IIl2(A) 


with I = {h^h), 1^1 = + ^2 and Let /x > 0 be a non-integer and 0 < s < 1, n 

is a non-negative integer such that n < < n + 1. We introduce the interpolation norms 

/ f°° _ \ ^ \ 

■= := , (4.1) 

where 

K(v,t) ;= ^^inf (||a - , 

Relevant Sobolev spaces are 

H>^{A) := {v G L^(A); < oo}> ■= {v e ^^(A); < oo}. (4.2) 

For /X > 0, it is known that .R^(R^) coincides with Ff^(M^). The following norms relation is 
useful in our analysis. 

Lemma 4.1. f77l[^7| / Let 0 < /x < 1 with fi ^ 1/2 and A he a domain in with Lipschitz 
boundary. Then H^{A) coincides with Hq{K) with equivalent norms. 

In the following, we shall give some useful results. 

Lemma 4.2. If dist{Ll'-,Llj) > IH for integer I >1, we have 

{D0 V,D^^^w) < j^^\\v\\a\\w\\a, veVi,weVj. (4.3) 



Figure 1: 
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Proof. Without loss of generality, we prove that (14.31) holds under the situation as the figure [T] 
shows. Cartesian coordinate x'Oy' is obtained by rotating xOy 9 angle counterclockwise. We 
set that Pi € O, and pj € such that the length of line segment piPj is equal to dist(n',np, 
the graphs of functions x' = Tii{y') and x' = Ti 2 {y') are parts of boundary Sfl' such that = 
{{x\y'):Ta{y')<x' < the graphs of functions Tji{y') and Tj 2 {y') are parts of boundary 

such that O'- = : Tji{y') < x' < Tj 2 {y')}, and the graph of x' = Toiu') is the 

perpendicular bisector of PiPj. Then we have Dg = _ooT)“,, = x'Dfo-, where -ooD’f, and 

x'D^ are the left and right Riemann-Liouville fractional derivative operator (see Appendix). We 
in the following only prove (14.31) for the case that there exists a ray in (cos 9, sin0) direction going 
through both and fl' from 14' to 14'-. Indeed, when the other case is true, {DgV, Dg_^^w) = 0 
and ()4.3p naturally holds. 

Denote 


M„/ = max y 

(x",y') e n' 

Mx' = max x‘ 

, y') e n' 


rriyi = min y'. 

(i'. y') e n', 

,y') 6 

nix' = aiin x'. 

(a:', y') e n' 


The shaded area is 


A 


<x'< Tj^2{y'),rny 


<y'< My'}. 


Let Aj = {{x',y')\x' < ro(2/')} n A and Aj = {{x',y')\x' > ro(y')} n A. Then we have 


{D^v, = (D^v, = mv, + mv, 

< \\D0v\\a^\\D^j^^w\\a^ + || D ^ n || A ,|| D ^+^ u )|| A ,. 

For {x,y) G Aj whose coordinate under x'Oy' is {x',y'), noting that supp(n) C 14', we have 


(4.4) 


d 


Dev{x,y) = -oc>D}f,v{x',y') =-^-ooDJ^ ,y') 


r(l — a) dx' 

1 d ^ 


< 


(x'— s) °‘v{s,y')ds 

p/. W'/ (x'- s)““n(s,?/')ds 

r(l - a) dx' Jri^iy') 

— -- / (x'- s)"(^+"^n(s,y')ds 

r(-«) JFniy') 

(^^)_(i+„) / \y(s,y')\ds 

JFu (v') 




rri2{y') 

'rn{y') 


1/2 


v^{s,y')ds 


(4.5) 


where in the fifth equality we have used the relation r(l — a) = —ar(—a), the first inequality is 
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by {x' — s) > lH/2 when s < Ti 2 {y') and the last by the Canchy-Schwarz inequality. Then 


pM , r^i2{y ) 

mv\\l < / dx' / dy' / V^{s,y')ds 

Jm^i Jniy, JTniy') 

, N f^y' r^a(y') 

< / dy' / v\s,y')ds 

dm, dViiiy') 


< H [[ v'^dx'dy' = 

J Jn' 


- 2 a ||,,||2 


(4.6) 


where the second inequality is by dx' < IH. Denote di as the diameter of D', and define a 
function in x' as 

0 , if x' > di] 


Hi{x') = { a,'(a-l) 

I 


if 0 < x' < di- 


For {x,y) G D' (whose coordinate is {x',y') under x'Oy'), by (IA.4D . we have 

v{x, y) = D""D^n(x, y) = -ooF>“" _ooD",?;(x', y') = Hi* _ooD'^,v{-, y'), (4.7) 

where v*^: denote the convolution product (see e.g., m)- Then by the Young Theorem (Theorem 
4.30 in [U), 

||<,y')ll[r.i,rd < \H;y')\k < ||77.|Ui(M)||_ooI?“,<,y')ll« < H^\\-ooD^^vi;y')\k. (4.8) 

Furthermore we have 


To' = 


< 


f/dx'dy'= f " M;y')\\fr,,iy^),r,,iy')]dy' 

n' J m^j 




TT2a\\ \\2 ^ TT2a\\ \\2 


\H°‘ (O) 


(4.9) 


where the second inequality is by Lemma 12.71 the last equality is by (j3.2p . Combining (|4.9I) with 
(SiSl), we obtain 

>IU- (4.10) 


mv\K ^ 


y ^ /1/2+c 
1 


Similarly we may also obtain 

WDe+MUi ^ (4-11) 

By Lemma 12.71 

\\D0 +^w\\a^ < ||a;||H-(o) ~ IklU, \\D^v\\Ai ^ \\v\\H<-{n) ~ ll^’IU- (4.12) 

Combining with ()4.4p , (|4.10p and ()4.1ip , we may obtain ()4.3p . □ 

Lemma 4.3. Let S C {1, 2,... , J} denote an index set such that D' n D'- = 0 for any i,j G S, 
i ^ j. Then for v = Yli&s Vi & Vi, we have 


{v,v)a < '^{vi,Vi)A- 

ieS 


(4.13) 
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Proof. Note that 


{v, v)a = B{v, v) = - {Dgv, Dg^^v)M{9)d9 + c(u, v), 

Jo 

it is easy to see that the lemma follows after we prove that 

< CY,ivi.Vi)A (4.14) 

i£S 

holds for any 9. So next, we give a proof of (j4.14p . It is easy to see that 

\{D^v,D^^^v)\ < Y. (4.15) 

k,j^S 

In fact, we may prove 

< cYiyi:yi)A (4.16) 

kjGS i^S 

through using the following two kind of inequalities: 

{D^Vk,D^^^Vj) < \\D^Vk\\\\D^+^Vj\\ 

< C\\vk\\A\\vj\\A<C\\vk\\A + C\\vj\\\, if < H (4.17) 

(by Lemma [22] and (|3.2p l: 

{D^Vk,D^^^Vj) < |o^||^^fe|UII^;ilU < j^:;^i\\vk\\A + \\vj\\A), 'li dist{n’^,n'j) >IH (4.18) 

(by Lemma l4.2p . 

Fixing i ^ S, we write 

Si = {j eS: dist^n'i, n'j) < H}, S'i = {j eS: dist{n[, n'j) > H}. (4.19) 

It is easy to see that 


the sum of the terms containing Vi on the left hand of (j4.16p 

< mv „+ E 

j&S jeS 

= Y ^ i(zi^„zi,v^,)i 

iGSi JgSi 

+ mvuD^^^v,)\ + mv„D^^^v,)\. (4.20) 

165' ies' 

We know that 

card{Si) < C, (4.21) 

where card{S) denotes the number of elements contained in 5, C is a positive constant dependent 
of ,1^. By (I4.2ip and (I4.17p . we have 

\{DgVj, D 0 _^_.jj.Vi)\ < Cllujll^ + terms which do not contain Uj. (4.22) 

j&Si j&Si 
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Figure 2: 


As figure [ 2 ] shows, L,, containing denotes the shaded bar area whose edges are parallel to 
the vector (cos(0), sin(0)) and go through Q'. It is not hard to see that {DgVi, Dg^^Vj) = 0 if 
Q'j Pi Lj = 0 , and for positive integer I, the number of the elements of the set 5- ^ 12 ^ fl Lj / 

0, IH < dist(f2'^, 12') < {l + l)H} is bounded by a positive C which is a positive constant dependent 
of And then using (I4.18p . we have 


< 

< 

< 


265 ' j€Sl 

Y mvuDi^^v,)\+Y E mvj,D^+.v^)\ 

i>i jes' i>i j€S' 


C 


5 _)_^ Ibill^ + terms which do not contain Vi 


1=1 


C||ui||^ + terms which do not contain Uj. 


(4.23) 


From above analysis, using (I4.17p and (|4.18D to enlarge the correspondent terms of the left hand 
of ()4.16l) . we may show that the inequality (|4.16l) holds with constant C being independent of 

e. □ 

Lemma 4.4. There exist linear operators Ki (i = 0,1,J), kq maps into Vq, Ki{i = 

1 ,..., J) maps and Vh into Hq{Q.[) and V) respectively, such that for any v G Hq{Q), 

j 

u = (4.24) 

i=0 


and 


J 


Y 11 ^*^ 

i=0 




J 


E11'^*'''11^1(0) 

i=0 


< c 



2 

Aqn)’ 


where the constant C depends on . 


(4.25) 


(4.26) 


11 















Proof. We employ the symbols used in the Lemma 3.12 of [28]. First let : Hq{Q,) —>■ Vh be a 
quasi-interpolation operator, which is defined as, for the node y in the coarse triangulation F//, 


ii^w)iy) 


0, y G dQ, 

\^y\~^ loj w{x)dx, otherwise. 


where uiy denotes the union of elements in F^f that share y, \uy\ denotes the measure of Uy. Let 

be the piecewise linear interpolation on the fine triangulation F/j, be the projector from 

into Vh, and 9i {i = 1,..., J) he the modified unit partitions of (see (3.7) in [28]) 

which are piecewise linear functions on Vh, satisfying ||V6*i(x)||oo < § and 

1<2< J 

supp( 6 'i) = ^i- ^ 

Take kq = and define /tj, i = 1, 2,..., J as follows: for w G Hq{Q), 

KiW := - Ko)w) + 9i{I - Q^)w. 


It is not hard to see that kq maps Hq{Q) into Vq, Ki{i = 1,...,J) maps Hq{^) and 14 into 
.^g(n') and Vi respectively, and (|4.24l) holds. Furthermore, define 4 := I^{9iQ^{I — KQ)v),Vi := 
9i{I — Q^)v, i = 1,2,... , J. It is easy to see that 


'^\\Kiv\f < ||kou|P + C'^|| 4||^ + C'^||4|P, 


i=0 


El 

2 = 0 


2=1 

J 


2 = 1 






2=1 


2 = 1 


| 2 _ 


In order to show that (|4.25l) and (I4.26p hold, it suffices to show that 


(4.27) 


(4.28) 


||kou||2 < C\\v\\^, 


(4.29) 


^ II4IP < C\\v\f,'^\\vi\f < C\\v\f, 
2=1 

Il^o^lli7i{o) - ^Il^lli7i(0)’ 


2=1 


and 


J 

El 


H 


J 


§ y ii'^iij7i{o)’ X/ 11^*11^1(0) - 
2 = 1 


< C 1 1 — ) ||u||£>l,o^, > . ll4ll£>Wo^ < C* I 1 + — ) 14 


H 


5 / ii"ii//i(o)- 


\w - Q’^wW + < Ch\\w\\^i^^y w G 


2=1 

It is known that 


By Lemma 3.6 in [25], we know 

lire - < CH\\w\\fji(^^y w G 

By a similar proof of Lemma 3.9 in [28| . we may get 

II^^''^IIj7i(o) ^ C'll^lli7i(0)> < C'lkll, 


(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 

(4.35) 
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here u; is a continuous piecewise quadratic function on F/j. 

Then ()4.29|) and (I4.31ji are a direct consequence of (I4.35p . meanwhile (|4.30p and (I4.32p can be 
proved by using (14.331) . (|4.34p and (j4.35p (Please refer to Lemma 3.12 in [28] for details). □ 
Next, we shall present the main results for the preconditioner system (|3.5p . 

Lemma 4.5. The operator is symmetric positive definite with respect to the inner product 
{■r)A. 


Proof. For the proof, we refer to section 2.3 in [28] . 

Lemma 4.6. For v G Vh, we have 

{PhV,v)A < C{v,v)a- 

Proof. By the definition of Pq, we have 

{Pqv, v)a = {Pov, Pqv)a < {v, v)a- (4.36) 

We recall the index set fi {i = 1,2,..., J^): for any j, k & fi with j 7^ k, 12) n 12'^, = 0, so it is easy 
to see that 

{PjV,w)A = {PiV,w)A, eVj,j e Ii,v eVh- (4.37) 

For any w € Vi, decompose it as tc = '^j ^ ^j- Then by ()4.37p . we have 

'^{PjV,Wj)A = '^{PiV,Wj)A = {PiV,w)A, v€Vh. (4.38) 

je/i jeii 

Taking Wj = PjV in the above equation, and using Cauchy-Schwarz inequality and Lemma 14.31 
we may obtain 

'^{PjV,PjV)A = {PiV,'^PjV)A 


je/i 


j&ii 


1/2 


< WPMUW'P^ Pjv\\A < \\Piv\\A [ 'P^\\Pjv\\A 


j&Ii 




which is 

It is easy to see that 
and 


^iPjV,PjV)A<\\Pi 

j^h 


Ull^. 


{PiV,v)A = {PiV,PiV)A < {v,v)a 


J 2 {PjV,PjV)A < \\PiV\\A < {v,v)a. 

iG/i 

Combining (14.421) with (|4.36p . we arrive at 

{Phv, v)a = {Pqv, Pqv)a + {Piv, Piv)a H-h {Pjv, Pjv)a 

< [Pqv, Pov)a + ^ (TjU, Pjv)A + • • • + 

j&h 

< {1 + J^){v,v)a. 

Then the lemma is proved. □ 


(4.39) 

(4.40) 

(4.41) 

(4.42) 


(4.43) 
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Lemma 4.7. For v G Vh, there exists a decomposition v = ^ such that 

j 


'^^{viTVi)A < C* ( 1 + 


i=0 


H 




(4.44) 


Proof. Let Ki be the operators as defined in Lemma 14.41 we now prove the decomposition v = 
Yli=o'^i with Vi = KiV may satisfy (I4.44p . Since ||tc||^ « ll'^llfl'“(n) ^ ^ suffices to show 
that 

J 


El 

i=0 


mv 




<C( 1 + ^ 1 \\v\\U 




(4.45) 


Actually, using the definition of the norms || • || • ||j^c.(Qp Lemma HTTl and Lemma we 

have 


Moreover, 


,7 J 

2=0 2=0 

^ /*oo 

y~] / 

«)iG77i(R2) 




ACyT; — Wj 


I .211 112_ 




< 


E 


inf 


J coo 


KiU — Wj 


211 112 \ 


-2«-l 


dt 




-2a-l 


j^O-^0 «)iG77l(n) 




dt 


-2a-l 


dt 


Jo «;€77l{f7) ^ V ^ ^ 

+ f) r J 4 n) ^ 


< 1 


( 1+ ^ j \M\H-‘in) 


1+ ^ j ll^liH“(n)- 


-2o-l 


l77“(0) 


< 


< 


r (I 

Jo w€H^(R^) ^ 


dt 


-2a-l 


dt 


V — W 


+ i'lkll|i(K2) 


t-^^-^dt 




I77“(R2) ~ IPll77“(0)- 

Combining (14.4711 with (14.4611 . we may get (14.4511 . □ 


(4.46) 


(4.47) 
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Lemma 4.8. For v G Vh, we have 


{v,v)a<C [1 + ^] {PhV,v)A- 


Proof. By Lemma W7f\ we know that there exists a decomposition v = ^ such that 

(j4.44p holds. Then by the Cauchy-Schwarz inequality and (14.441) . we have 

J J 

(v,v)a = y^^(v,vi)A = 'y](Piv,vi)A 

i=0 i=0 

/ J ^ \ 

'^[PiV, PiV)A 


< 


^*=0 
^ J 


^j=0 


\ V 2 . J s 1/2 


vi=0 


vi=0 


< {PhV, vfl^ (^C ^1 + y ^ (u, 


1/2 


which is Lemma 14.81 □. 

Theorem 4.9. The condition number of P^ satisfies 

cond{Ph) < C n + y ) . 


(4.48) 


Proof. The theorem follows from Lemma 14.61 and Lemma 14.81 □ 


5 Numerical results 

We test PCG method with as a preconditioner. The tests are carried out by using Matlab 
software. The stopping criterion is 

\\u^ — tt^“^||oo < 10“®, 

where is obtained from by one step PCG iteration, and = 0. Here we only test the 
case that H is a square domain. We take the triangulations Th,^^ as uniform ones as in figure [3] 
shows. With viewing the one in figure [3] as the coarse triangulation Th, the subdomains and 
are taken as the darker and lighter shaded squares respectively. 

We give two examples in this paper: one is with the probability measure M having a discrete 
form and the other is with M being a continuous function. Table [1] and Table [2] list the iterative 
steps of our PCG methods for Example 15.11 and Example 15.21 respectively. It can be seen that 
when we fix ^ = constant, for instance -f- = 8, the convergence rate of our PCG method is 
optimal, i.e., the iterative steps of the PGG method almost keep to be a constant, which coincides 
with our theoretical results in this paper. Prom the tables 1-2, we also find that small overlapping 
case is a good choice for our domain decomposition methods for the SFPDEs, which is same as 
the overlapping domain decomposition methods for the integer order PDEs. 
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n, 

ti' 

n 


Figure 3: The trangulation and the subdomains 


(h,H) 

II 

5 = 2h 

II 

(A. 'A) 

13 

13 

14 

(tM’ Te) 

12 

11 

12 

(2 'i 

V 256 > 82 ) 

10 

10 

11 


Table 1: The iterative steps of PCG method for Example 15.11 

Example 5.1. Let H = [0,2] x [0,2], the equation is 

-+ yDl,^)u = fix, y), (5.1) 

where fix,y) is taken such that the solution is u = (( 2 x — x^)( 2 y — y^))^. 

Example 5.2. Let Q = [0,2] x [0,2], the equation is 

- = fix,y), (5.2) 

where M(0) = 1, 0 € [0, 27r), and fix,y) is taken as in Example 15.11 



Appendix 

In this appendix, we shall present the definitions of the fractional integrals and derivatives in one 
dimension case, and their corresponding properties, which are used in this paper. The related 
results can also be found in [2Tl[2l]. For function u(x),x € (a,&), where a may be taken as —oo, 
and b as +oo, the left and right Riemann-Liouville fractional Integrals of order y > 0 are defined 
respectively by 

aDf^vix) = J xD~^vix) = J xV~^vis)ds. (A.l) 


(h,H) 

II 

6 = 2h 

II 

(A. 'A) 

13 

13 

13 

(m’ Te) 

13 

13 

13 

( 2 2 \ 

V 256 ’ 82 > 

13 

12 

13 


Table 2: The iterative steps of PCG method for Example 15.21 
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The left and right Riemann-Liouville fractional derivatives of order /r > 0 are defined respectively 



where n is an integer such that n—l<^<n. For p > 0, A: is an integer satisfying k — 1 < p < k. 
If aDxv(x) is integrable, then 





(A.4) 


(jA.4p is from (2.108) in [2T]. 
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